Low-energy quenching of positronium by helium 



> 

o 

m 

(N 



C/3 

1j 



Sadhan K. Adhikari, P. K. Biswas, and R. A. Sultanov 

Instituto de Fisica Teorica, Universidade Estadual Paulista, 01.405-900 Sao Paulo, Sao Paulo, Brazil 

(February 2, 2008) 

Very low-energy scattering of orthopositronium by helium has been investigated for simultane- 
ous study of elastic cross section and pick-off quenching rate using a model exchange potential. 
The present calculational scheme, while agrees with the measured cross section of Skalsey et 
al, reproduces successfully the parameter ^Z^^, the effective number of electrons per atom in a 
singlet state relative to the positron. Together with the fact that this model potential also leads 
to an agreement with measured medium energy cross sections of this system, this study seems 
to resolve the long-standing discrepancy at low energies among different theoretical calculations 
and experimental measurements. 
PACS Number(s): 34.10.+X, 36.10.Dr 
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Studies on positronium- (Ps-)impact scattering 
have gained momentum these days due to the avail- 
ability of ortho-Ps beam in the laboratory and its vast 
apphcational potential coupled with the present in- 
adequate and inconclusive understanding of its inter- 
action dynamics with matter jl]. Ps scattering by 
neutral targets has posed a challenge to theoreticians 
on a proper accounting of experimental data as most 
existing theoretical works disagree with the major ex- 
perimental trend. 

The discrepancy figures out prominently in the Ps- 
He system where there are many theoretical and ex- 
perimental studies. The medium-energy experimental 
cross section shows a declining trend with decreasing 
energy Q from a peak around 20 eV for Ps-He scat- 
tering. Similar trend is also observed in PS-H2 and Ps- 
Ar systems [^,||. This trend, which is supported by 
the recent measurement of Skalsey et al [|| , could not 
be reproduced in most theoretical predictions 
At very low energies, these theories and experiments 
[p|-^, p^pl| ] on the Ps-He system are inconsistent with 
each other and also among themselves. For illustra- 
tion, the zero-energy cross sections calculated on Ps- 
He by different authors vary from 3.3A^ to 16.54 
11^ while the measured values range from 2.3±0.4 
A2 p 0.915 eV) to 11 ± 3 A2 (between to 0.3 
eV) Q . Pointing out the very reactive nature of Ps 
scattering and its associated convergence difficulties, 
a prescription for the generation of nonlocal model 
exchange potential has been advocated recently and 
applied successfully to different electron-impact (tar- 
gets: H, He) and Ps- impact (targets: H |l^Jl^, He 
[0, H2 III], Ar, Ne (ID) scatterings problems using 
static-exchange to three-Ps-state models. The three- 
Ps-state calculation for Ps-He predict a further lower 
zero-energy cross section of 2.42 A^ |Q. 

In this work we shed light on the abovementioned 
discrepancy in Ps-He system in conjunction with a 
determination of the parameter ^Z^^ which denotes 
the effective number of electrons per atom in a singlet 



state relative to the positron. The study of this pa- 
rameter is supposed to provide a stringent test for the 
model potential. The incident ortho Ps(ls) atom in 
a triplet state with a lifetime of 142ns can decay into 
three photons and is more stable than its para coun- 
terpart in a singlet state with a lifetime of 0.125ns for 
a two-photon decay mode. However, in its interaction 
with matter, the positron of Ps can find an atomic 
electron in a spin-singlet state and the two can be 
annihilated by a two-photon decay mode without re- 
ally forming a para Ps atom by electron exchange. 
This process is termed pick-off quenching. From the 
experimental pick-off quenching rate the parameter 



^eff 



can be extracted. 



Theoretically ^-^gff can be 



calculated from the wave function of the Ps-He sys- 
tem ^(rp, Sp] ri, Si; r2, S2; ra, S3) where r and s refer 
to position and spin, the suffix p refers to the positron 
and « = 1, 2, 3 refer to the electrons. Following Barker 
and Bransden |^ 13|, the amplitude for finding the 
positron and one of the atomic electrons in a relative 
singlet state is 

*(i-p; ri; r2,S2; rs, S3) =< Xo(sp, si)|* >, (0.1) 
where xo is the singlet wave function. The parameter 



^Zg5 = 3^ y drpdridr2dr3(5(rp - ri)|$| 



(0.2) 



spin 



The factor 3 appears as each of the three electrons of 
the Ps-He system contributes equally to ^^gff- 

Unlike the scattering cross sections, which are de- 
termined from the asymptotic part of the Ps-He wave 
function, the parameter ^^gff is sensitive to the Ps-He 
wave function at short distances and its correct eval- 
uation in a theoretical calculation should provide a 
sensitive test about its realistic nature. There is con- 
siderable discrepancy between theory and experiment 
in the value of the parameter ^ Z^^. The experimen- 
tal measurements have yielded ^-^gff = 0.108 ± 0.01 
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0.135 ± 0.068 and ^Z^q = 

0.25 ± 25% ]7y2l[|, whereas different static-exchange 
calculations have yielded values ranging from 0.02 to 
0.1 |||,0,|l^. Compared to other exchange potentials, 
the present scheme leads to substantially weaker re- 
pulsive exchange potential and it is expected that the 
present scheme will lead to a larger value of ^Z^ff in 



Ps-He system as is demanded by experiments 19 21 



In the static-exchange approximation the Ps-He 
wave function is represented by the following an- 
tisymmetrized product of the internal wave func- 
tions of Ps(ls), (/)ps(ri — Tp), and singlet He(lsls), 
0He(r2,r3), with a wave function of relative motion, 
Fk(R), and a suitable spin function: 

), (0.3) 

where R = (ri + rp)/2, t = ri — rp, A is the anti- 
symmetrizer, k the incident Ps momentum, and x the 
spin function. 

On expanding -F(R) in partial waves 

oc 

Fk(R) = ^(2L + l)(fci?)-iFi(i?)PL(cos0), (0.4) 

where 9 is the angle between k and R, the follow- 
ing integro-differential equation is obtained from the 
Schrodinger equation: 



( — 



L{L + l) 



FLiR) 



i?2 

T4(i?, R')FL{R')R'^dR', 



(0.5) 



where we use Vl{R,R') is the nonlocal exchange 
potential. The asymptotic boundary conditions for 
Fk(R) and Fl{R) are given by 



Fk(R) 

Fl{R) 



^R^oo exp{tkRcos9) + f{e f'^^^]^^\ (0.6) 



R 



sin(fci? - i7r/2 + (5l), 



(0.7) 



where (5l is the scattering phase shift, and the scat- 
tering amplitude f{9) is given by 

m-±{2L^lf^^^^^4^Pd^o.9). (0.8) 



L=0 



k 



The total elastic and momentum transfer cross sec- 
tions are given by 



<7el {k^ 
cr„(fc^ 



\I{0)?dn, (0.9) 
|/(6i)|2(l-cos6')df], (0.10) 



respectively. Some of the experiments provide 
only low-energy momentum-transfer cross section and 
hence we also calculate this observable in this study. 



We employ He(lsls) wave function of the following 
form 



^'ffe(r2,r3) = U2(r2)w3(r3) 
Mi(r) = y^aB»0Kt(r), 



(0.11) 
(0.12) 



with 4>i^i{Y) — exp(— aKi'')^oo- In the present calcula- 
tion we use an accurate two-term parametrization of 
( 3.12 ) [^ . In momentum space the model exchange 
potential has the following form jl^ : 



B(kf,ki) 



EE 



3 KK' 



J exp(iQ.r)<?!)«,j(r)dr 

J (/)J,,(t)exp(iQ.t/2),^p,(t)dt, (0.13) 



with Q = k; - kf , V{p, q) = -B{p, q)/(27r2) and 



(0.14) 



where 0Kj"(r) is the Kth function of the jth elec- 
tron for the atomic ground state, and (?!)ps(t) = 
P^^^ exp(— /3i) / The direct potential for this prob- 
lem is zero, and there is a change of sign in the spin- 
triplet Ps-He potential below. The partial-wa ve c on- 
figuration space nonlocal potential of Eq. (|0. 
given by 



Vl{R,R') 



00 />oo 



p^dpq^dq 

J Jo Jo 
X jL{pR)VL{p,q)jL{qR') 



(0.15) 



and 



VL{p,q)^-TT^ J dxPL{x)V {p,q), (0.16) 

where x is the angle betwee n p a nd q. Although the 
model exchange potential ( |0.16 ) has been found to 
be satisfactory for calculating scattering cross sections 
[|l4| , it is interesting to investigate if it is also effective 
in the calculation of finer scattering observables, such 
as the parameter ^Z^q. A scheme for the calculation 
of ^Z^Q for He wave functions of type ( |0.11 ) is given 
by Fraser and Kraidy [p^ and we employ the same in 
the present calculation. 

First we show our results for the low-energy elastic 
and momentum-transfer cross sections in Fig. 1 to- 
gether those obtained from other theories and experi- 
ments. The triplet Ps-He scattering length in this case 
is 0.87ao, compared to 1.39ao obtained by Drachman 
and Houston ||]. The discrepancy among various re- 
sults is apparent in this plot. The three experimental 
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results for cross sections shown by solid circles 

vary from 2.3 [|, through 7.45 to 11 A^ 

at 0.9 eV, eV, and 0.15 eV, respectively It is 
difficult to reconcile these three experimental results 
in a theoretical model. Previous static-exchange cal- 
culations |0-||,|l^ except those of Rcf. [|l^ all tend 
to support the largest cross section of Ref. ||l^. The 
model calculation of Drachman and Houston |^ de- 
noted by a cross in Fig. 1 is consistent with the exper- 
iment of Ref. . This discrepancy has been partially 
resolved in the study of Ref. , where it has been 
demonstrated that the present model exchange poten- 
tial is unique in being able to reproduce experimental 
cross sections j2]j^ upto medium energies (about 60 
eV) fairly well. Other theoretical calculations are un- 
able to reproduce ||^,^,|l^ the experimental trend of 
total cross section at different energies with a mini- 
mum around 5.1 eV. The present elastic (full line) and 
momentum-transfer (dashed line) static-exchange cal- 
culations are consistent with the experiment of Skalsey 
et al 01. 



Figure 1 
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Fig. 



1. Angle-integrated Ps-He cross section at 
low positronium energies: present momentum trans- 
fer from static-exchange model (dashed line); present 
elastic from static-exchange model (full line); elastic 
from static-exchange of Ref. j?) (dashed-dotted line) ; 
of Ref. |l|l (dashed-double-dotted line); of Ref. B 
(dashed-triple-dotted line) ; theory of Ref. |^ (cross) ; 
experiments at eV, 0.15 eV, and 0.9 eV of Refs. 
|pT| , P^ j4|], respectively (solid circle). 

Next we perform an S-wave calculation for the pa- 
rameter ^Z^Q. In Fig. 2 we plot the results for the 
present calculation of ^Z^q at different energies. For 
comparison we also plot the results of previous calcu- 
lation by Barker and Bransden 0| and by Fraser and 
Kraidy p3[ | and the existing three experimental data. 
Although Drachman and Houston g did not calcu- 
late ^^gff at different energies, their low-energy value 
of 0.1 is in reasonable agreement with the present re- 
sult of 0.11 and experiment of Refs. The three 
experimental results with error bars cover the range 



of 0.07 to 0.31 for ^Z^q. Of the three experiments the 
one by Duff and Heymann with the smallest error 
bar might be the most accurate. 

The much too small values of ^^gff obtained in 
previous calculations [^jl3| seem to be a consequence 
of a much stronger (exchange) repulsion in the elas- 
tic channel of these models. This is reflected in the 
zero-energy cross section or the scattering length of 
these calculations. For a repulsive potential the low- 
energy cross section increases with repulsion, conse- 
quently, the previous calculations have led to unusu- 
ally large triplet scattering lengths compared to the 
present work. This is most clearly exhibited in a cor- 
relation exhibited in Fig. 3 where we plot ^Z^q versus 
triplet scattering length of different calculations. The 
larger the scattering length the smaller is the ^Z^g. 
This correlation is similar to different correlations ob- 
served in the study of Ps-H scattering in Ref. ||l^ . 

Figure 2 
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Fig. 2. The parameter ^Z^q at different positron- 
ium energies: calculation including angular momenta 
L = 0, 1, 2 of Rcf. |l| (dashed-dotted hue), of Ref. § 
(dashed line); calculation for L = of present model 
(full line); the experimental points denoted solid cir- 
cle, diamond, cross taken from Refs. p9|-pH|, respec- 
tively. 

We next comment on two aspects of the present 
calculation. First, we used a two-term helium wave 
function. We also repeated our calculation with the 
one-term helium wave function of Ref. ||^ and the five- 
term wave function used in Ref. [Q. The results for 
both the cross section and ^Z^q suffer insignificant 
change with the change of wave function. For ^ Z^q, 
the different results are within the error bar of Ref. 
Jigf ; for cross section they are also within the error 
bar of Ref. jj]. Hence we do not believe the present 
results to be so peculiar as to be of no general validity. 
Secondly, we performed a L = calculation for ^Z^q. 
At the experimental energies less than 0.03 eV, the 
effect of higher partial waves is practically zero (well 
within the error bar of Ref. fl^ ) ; at 1 eV this effect is 
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quite small. 



Figure 3 
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Fig. 3. The parameter ^Z^g versus scattering 
length of various calculations denoted by solid circle 
and labeled by reference number j| . The present 
result is labeled by [PR], full line denotes a linear fit. 

From the consideration above, we believe the 
present static-exchange calculation as well as the pre- 
vious calculations using the same exchange po- 
tential provide a realistic account of very low-energy 
Ps-He scattering. However, the precise agreement of 
the static-exchange calculation with experiment is ex- 
pected to be incidental. For a complete understanding 
of this problem higher excited states of both Ps and 
He should be incorporated in the model. The inclu- 
sion of Ps excitation channels has been found |Q to 
decrease the low-energy cross sections and we might 
need to refit the low-energy cross sections by chang- 
ing the parameters a and/or f3 of the potential in Eq. 
( |0!l4| ), as in Refs. fl|,|l|. 

In conclusion, we have used a recently suggested 
nonlocal model exchange potential [ p^p^[p^Jl7[ | and 
applied it to the study of Ps-He scattering at low ener- 
gies. We have critically examined the static-exchange 
calculation to see if it can account for satisfactorily 
[0 the measured cross sections of Refs. and 
the 'measured' value of the parameter ^Z^q of Refs. 
^ , pO| . The present calculation is in reasonable agree- 
ment with the calculation of Ref. However, it is 
difficult to reconcile the present calculation with the 
experiment of Ref. |2^] . In a previous study , the 
present exchange potential has been found to repro- 
duce the low- to medium-energy cross sections of Refs. 

well. This coupled with the present study seems 
to indicate that the ^Z^q measurement of Ref. 
and the high- and low-energy cross section measure- 
ments of Refs. dl and jj] are consistent among each 
other as well as with the present calculation which 
possibly provides a faithful description of low-energy 
Ps-He scattering. We observe a correlation between 
^ZgQ- and triplet scattering length of various calcula- 
tions (Fig. 3), which demonstrates that the smaller 
the scattering length the larger is the value of ^Z^g. 
This correlation is similar to different correlations ob- 



served between the low-energy Ps-H observables re- 
cently @. 
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